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Abstract 

We develop a coupled-mode theory for spatial gap solitons in the one-dimensional photonic 
lattices induced by interfering optical beams in a nonlinear photorefractive crystal. We derive a 
novel system of coupled-mode equations for two counter-propagating probe waves, and find its 
analytical solutions for stationary gap solitons. We also predict the existence of moving (or tilted) 
gap solitons and study numerically soliton collisions. 
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It is well known that a periodic modulation of the optical refractive index not only 
modifies the spectrum of linear waves, but also strongly affects nonlinear propagation and 
self-trapping of light . Recently, formation of spatial solitons in reconfigurable photonic 
lattices in photorefractive materials has been demonstrated experimentally in one- and 
two- [5( dimensional geometries. In this case, strong electro-optic anisotropy of a photore- 
fractive crystal was employed to create the lattice by interfering laser beams with ordinary 
polarization, while the solitons were observed in the orthogonally polarized mode. 

Periodically modulated nonlinear systems can also support self-trapped localized pulses 
or beams in the form of gap solitons (GSs), which are hosted by a gap of the system's linear 
spectrum, induced by the resonant Bragg coupling between the forward- and backward- 
propagating waves [&]. A notable property of the GSs are that they can form in both 
self-focusing and self-defocusing media. 

Travelling temporal-domain GSs have been observed experimentally in material Bragg 
gratings written in silica fibers Q|. The concept of a spatial-domain GS was proposed Js[ £fl 
and further elaborated in waveguide settings jlfl llll ll^. Experimentally, spatial GSs were 
demonstrated in waveguide arrays and optically- induced photonic lattices 14]. 

The simplest and most ubiquitous description of the GSs is provided by the coupled- 
mode theory (CMT), which amounts to the derivation of a system of coupled nonlinear 
propagation equations for the forward and backward waves p. In this Letter, we derive 
a novel type of the CMT model for spatial GSs in the optically-induced photonic lattice. 
We find exact analytical solutions for GSs in the framework of this model, and compare 
them with the results obtained numerically by solving the full nonlinear model. We thus 
identify a parameter region where the couple-mode approximation is accurate enough, and 
in that region we use the model to predict new features such the existence of the families 
of moving (or tilted) gap solitons and to study their interaction Actually, the new model 
derived in this work may have more general purport than just an asymptotic approximation 
for the nonlinear photonic lattice, as it is the first coupled-mode system that accounts for a 
saturable nonlinearity. 

Following Refs. P, Q], we consider the evolution of an extraordinarily polarized probe 
beam propagating through a periodic structure written by counter-propagating ordinary- 
polarized plane-wave beams in the photorefractive medium. As mentioned above, the electro- 
optic coefficients in the crystal are substantially different for the orthogonal polarizations, 
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therefore the grating, created by interference of the counter-propagating beams, is essentially 
one-dimensional (uniform along the propagation axis z). The interference of the plane 
waves with a wavelength A creates an intensity distribution, I g (x) = IqCos 2 (Kx), with 
K = 27rn A~ 1 sin 9, where 9 is the angle between the Poynting vectors of the plane waves 
and the z axis, and no is the refractive index along the ordinary axis. Provided that the 
intensity of the probe beam, \E\ 2 , is much weaker than that of the grating, Jo, one may 
neglect the feedback action of the probe beam on the grating (a perturbative calculation 
within the framework of an extended model, that includes equations for both the lattice- 
forming and probe fields, shows that the feedback exactly cancels out at the first order in 



1/Jo). Thent 



equation , whose normalized form is 

.dE Id^E E 

% dz + 2 dx 2 1 + J cos 2 (Kx) + \E\ 2 ~ U 

To derive the corresponding equations for the coupled forward and backward waves, we 
approximate solutions to Eq. (JTJ) by 

E(x, z) = u(x, z)e iKx + v{x, z)e~ iKx , (2) 

where u and v are slowly varying [in comparison with the carrier waves exp (HKx)] envelopes 
of the forward and backward modes. Substituting the expansion (J2J) into Eq. ((TJ), we perform 
the Fourier expansion with respect to exp (±iKx) and, in the spirit of the CMT approach, 
keep only the lowest-order harmonics. Eventually, this leads to the following CMT equations 
(one of which is linear): 

d d 
i — (u-v)+iK—(u + v) (3) 

oz ox 

2(U ~ V) 0, (4) 



re evolution of the local amplitude E(x, z) of the probe beam obeys the known 



i (i + 1« - tf ) + 1 + 2 (M 2 + M 2 ) 

i—(u + v)+iK—(u-v)=0. (5) 

Equations (jlj and (0) constitute a new CMT model with the saturable nonlinearity. It 
contains one irreducible parameter Jo, while K can be absorbed by rescaling of x. 

In the physically relevant case, the photonic-lattice intensity is large, i.e., J 3> 1, \u\ 2 , \v\ 2 , 



hence the square root in Eq. (jlj) may be approximated by J hi)- + \u — v | 2 ), except for a 
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vicinity of point (s) where w = u — v vanishes. Using this approximation, and eliminating 
(u + v) by means of Eq. (jSJ), we reduce Eq. (jlj) to an equation for the single function w, 

K^ + ^l\ ~ 1=0. (6) 




dz 2 

We have checked the accuracy of the simplified equation ©, comparing its analytical solu- 
tions for solitons (see below), and their stability, vs. direct numerical solutions of Eqs. (J3J 
and A conspicuous difference appears only in the region of Jo ^ 3, where the CMT does 
not provide for an adequate approximation anyway. 

Stationary solutions of Eq. (0) are sought for as w(x, z) = e lQZ W(x), where a real function 
W obeys the equation d 2 W/dx 2 = —dU c s/dW, with 

U«(W) = ^- 2 w> + (Vl+W - l) . (7) 

Further, Eq. (JHJ) shows that, for the stationary solutions, 

< . 1 



W(x) ± 



q dx 

If the propagation constant q belongs to the interval (which is, actually, the bandgap) 



(8) 



< -q < 2/yj = Q, (9) 

the potential (jJJ has two symmetric minima, giving rise to GS solutions that can be written 
in an implicit analytical form, 

( KW'f = - {qWf + U\q\/Jl^) (Vl + W 2 - l) . (10) 



As it follows from Eq. (fTU|) . the soliton's squared amplitude is W^ ax = (8/g 2 /o)(2 — 
y/To\q\). In the small-amplitude limit, i.e., for < e = 2(2 — |<z|v%) «C 1 [cf. Eq. 
©], the GS asymptotically coincides with the conventional nonlinear-Schrodinger soliton, 
W{x) = sech ^e\q\/2K~ 1 I 1 ^ 4 x^) . In the other limit, \q\ — > 0, the soliton assumes 

a "compacton" shape: \q\\/l\ ) W(x) = 4 cos 2 (\q\K~ 1 lQ 1 ^x s J , if |x| < ttKIq^/ (2\q\), and 
W(x) = otherwise. However, the conditions under which the CMT equations were derived 
above do not hold in the latter case. 

To verify the validity of the CMT approximation in the present setting, in Fig. ^ we 
compare the analytical soliton solutions based on Eq. (fTU|) and ones found numerically from 



Eq. (JIJ . The comparison is presented in terms of a global characteristic of the soliton family, 
viz., the integral power, N = f_™ \E(x)\ 2 dx, vs. the propagation constant q; an example 
of the comparison of the soliton's shapes is also included. Naturally, the approximation is 
appropriate sufficiently close to the gap's edge. We also note that the negative slope of the 
N(q) dependence suggests stability of the entire soliton family as per the Vakhitov-Kolokolov 
criterion, i.e., the absence of real eigenvalues in the spectrum of small perturbations around 
the soliton [15]. However, the solitons may be subject to instabilities with complex eigen- 
values. Detailed numerical analysis shows that, in the case of large Jo, when the the CMT 
approximation is relevant, a stable part of the soliton family in the bandgap (jHJ) is limited 
to a sub- gap, 

Q-AQ<-q<Q. (11) 

The width of the stability sub-gap, AQ rs 0.21, is nearly constant for Jo > 20 [for Jo = 25.5, 
which is the case shown in Fig. ^ the GSs are stable in « 53% of the interval (JSJ); we mention, 
for comparison, that in the standard GS model with the cubic nonlinearity, the stable 
part occupies ~ 51% of the bandgap is |3,[l2|]- Unstable solitons [which actually occupy 
that part of the bandgap (Q where the CMT approximation is irrelevant] are destroyed by 
perturbations. 

The CMT approximation opens many ways to investigate novel phenomena which may be 
difficult to study directly within the framework of the full model. An issue of obvious interest 
are moving (or tilted) GSs, of the form (u, v) = e iqz (U(x — cz), V(x — cz)). Our analysis 



shows that the bandgap Q for the tilted solitons shrinks to < —q < Q = Qyl — (c/K) 2 , 
and it does not exist for c 2 > K 2 . The reduced gap is completely filled with moving gap 
solitons which are stable in the sub-gap Q — AQ < —q < Q, cf. Eq. flllj) . with the width 
AQ that does not depend on the tilt c, up to the accuracy of numerical results (the latter 
fact resembles a known stability feature of the GSs in the conventional CMT approach }l7|). 

The stability of the tilted solitons suggests to consider collisions (intersections) between 
them. The most important characteristic of soliton interactions is elasticity. Our numerical 
simulations demonstrate that the interaction is quite elastic, unless the solitons are taken 
close to the instability border [which is — q = Q — AQ, see Eq. flllj) ]. In the latter case, there 
is a specific structure of regions of elastic and inelastic collisions in the plane (ci, — C2) of the 
tilts of two colliding solitons, as shown in Fig. El In the case of inelasticity, the colliding 
solitons suffer strong perturbations. 



In conclusion, we have developed the coupled-mode theory for spatial gap solitons in 
one-dimensional photonic lattices written optically in a photorefractive medium. We have 
derived a novel type of the coupled-mode equations and found their analytical solutions 
for gap solitons. In the case when the coupled-mode system correctly approximates the 
photonic lattice, we used the model to investigate new features of such systems, such as the 
existence and stability of moving gap solitons and their collisions. The latter results suggest 
new experiments for optically-induced lattices in photorefractive media. 

This work was partially supported by the Australian Research Council. B.A.M. appreci- 
ates hospitality of the Nonlinear Physics Center at the Australian National University. 
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FIG. 1: The integral power of the soliton vs. the propagation constant [in the first finite gap of 
Eq. ((H)] at Iq = 25.5 and K = 0.5 . The solid and dashed lines show, respectively, the numerical 
result, and the analytical one, as obtained from Eqs. (|lUj) . (JSJ), and ©. Examples of numerical (a) 
and analytical (b) profiles of a soliton taken close to the edge of the gap (q = 0.39, point A) are 
shown in the bottom part of the figure. 
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FIG. 2: Regions of elastic and inelastic head-on collisions between two solitons, with the tilts c\ > 
and C2 < 0, and \q\\ = \q2\ = 0.2 (close to the instability border, which is \q\ = Q — AQ ?a 0.19, in 
this case), at Iq = 25.5 and K = 0.5. 
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